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ABSTRACT 



In this Letter we establish clear evidence for the resonant absorption damp- 
ing mechanism by analyzing observational data from the novel Coronal Multi- 
Channel Polarimeter (CoMP). This instrument has established that in the solar 
corona there are ubiquitous propagating low amplitude (^1 km s~^) Alfvenic 
waves with a wide range of frequencies. Realistically interpreting these waves 
as the kink mode from magnetohydro dynamic (MHD) wave theory, they should 
exhibit a frequency dependent damping length due to resonant absorption, gov- 
erned by the TGV relation showing that transversal plasma inhomogeneity in 
coronal magnetic flux tubes causes them to act as natural low-pass filters. It 
is found that observed frequency dependence on damping length (up to about 
8 mHz) can be explained by the kink wave interpretation and furthermore, the 
spatially averaged equilibrium parameter describing the length scale of transverse 
plasma density inhomogeneity over a system of coronal loops is consistent with 
the range of values estimated from TRACE observations of standing kink modes. 

Subject headings: niagnetohydrodynamics (MHD) — plasmas — magnetic fields — 
Sun: corona — Sun: oscillations 
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Introduction 



Recently, there has been much interest generated by t he observation o f ubiq uitous 
propagating Alfvenic waves in the solar corona detected by 



Tomczyk et al 



( 120071 ) using 



the innovative Coronal Multi-Channel Polarimeter (CoMP) instrument. The Alfvenic 
properties of these waves are undeniable since they have a phase speed of about 1 Mm s~^, 
the velocity components are perpendicular to the direction of magnetic field lines and they 
exhibit only very weak intensity fluctuations, suggesting an almost incompressible wave 
mode. Detected in Doppler velocity, they have a low amplitude of approximately 1 km 
s~^ peak to peak corresponding to a loop displacement amplitude of only 48 km and peak 
power at a period of about 5 minutes. The small spatial extent of the displacement explains 
why these waves could not be detected previously by the Transition Region and Coronal 
Explorer (TRACE) since it has a resolution of about 730 km. However low amplitude, 
about 2 km s~^ peak to peak, Alfvenic waves in coronal loops, also with a peak power 
period of about 5 minutes, have been detected before in Doppler velocity by the E UV 



Imaging Spectrometer 



Van Doorsselaere et al. 



IS) i nstrument on Hinode, e.g.. 



Erdelyi fc Taroyan 



(120081) and 



( I2OO8I ). The consistent observed properties sug gest that the two 



i nstru ment s have detected the same type of propagating coronal wave. 



( 120081 ) and 



Van Doorsselaere et al. 



Erdelyi fc Taroyan 



( I2OO8I ) interpreted these trans versal coronal waves a s the 



kink mode from magnetohydrodynamic (MHD) wave theory (see 
a detailed discussion of this wave mode). 



Goossens et al. 



2009 



for 



2. Kink waves in coronal loops 



When t he first EUV post-flare transversal coronal loop os cillations were directly imaged 



by TRACE (JAschwanden et al. 



1999 



Nakariakov et al. 



19991), it was clear that the simple 



model of Alfven waves originally derived by Alfven (1942) based on an infinite homogeneous 
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plasma in planar geometry could not accurately describe the oscillatory properties of 
coronal loops, which are finite closed magnetic fiux tubes with inhomogeneous plasma 
structure. Coronal observations show the EUV intensity of oscillating loops is greater than 
the surrounding medium, indicating the internal plasma density is greater inside the loop 
than in the ambient plasma. Essentially, the simple homogeneous plasma model by Alfven 
would not permit us to attempt magnetoseis mology of the solar atmosphere to determine its 



fine structure (see e.g. 



Banerjee et al. 



20071 . for review) and would not allow us to estimate 



the contribution of wave heating to the coronal plasma, since the most likely mechanisms 
proposed thus far require plasma inhomogeneity in the direction transverse to the directio n 
of th e magnetic field, e.g., phase mixing and resonant absorption (see e.g.. 



Aschwanden 



2004 . for review of various proposed damping mechanisms). Due to the observed plasma 
inhomogeneity of coronal loops theorists concluded that the transversal oscillations detected 
with TRACE were best described as standing kink waves from MHD theory. This mode is a 
bulk transversal oscillation of magnetic fiux tubes that have different internal and external 
Alfven speed equilibria, primarily due to density variatio n in the direction transverse to 



the direction magnetic field (e.g.. 



Edwin &: Roberts 



19831 ). Since spectroscopic EUV line 



ratios can be imp le mented to determine coro nal d ensity using the CHIANTI dat abase (e.g. 



Landi et al. 



20061), lErdelvi fc TarovanI (120081 ) and 



Van Doorsselaere et al. 



(120081) exploited 



the spectral capabilities of EIS to determine that the coronal loops in which they detected 
low amplitude Alfvenic waves had greater internal plasma density than external, naturally 
leading them to the kink mode interpretation, in agreement with the previous consensus 
regarding transverse waves in coronal loops observed by TRACE. 



Notably, the observed coronal loop standing kink wa ves observed with 



Aschwanden et al. 



'RACE are 



20031) and the 



heavily attenuated in 1 — 4 fundamental periods (see e.g. 

most likely physical mechanism for this proposed thus far is resonant absorption caused 

by naturally occurring gradients in the transversal Alfven speed of coronal loops (see 
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e.g., iGoossena l2008l for review). Although resonant absorption is still a ro 



3ust damping 



mech anism in more complex geometries, e.g., multi-thread loop structures (JTerradas et al. 



20081 ). as an illustration of the basic process we assume a coronal loop to be a cylindrical 
magnetic flux tube. Then the kink phase speed in the long wavelength limit for an axially 
homogeneous magnetic cylinder with differing internal and external plasma density and 
constant magnetic field strength is 



where pi and pe are the internal and external plasma density, B magnetic field strength 
and p magnetic permeability. Equation ([T]) shows that the kink speed is intermediate 
between the external and internal Alfven speeds. Regarding kink waves in coronal loops, 
it is observed that pi > pe, so naturally occurring density gradients between the internal 
and external plasma will cause the local Alfven frequencies in this intermediate layer to 
match the global kink frequency at some magnetic surface, thus causing a resonance. This 
generates azimuthal Alfvenic motions in the resonant surface which grow in amplitude, 
damping the global kink mode as energy is transferred to this localized mode. For kink 
waves excited by a broadband disturbance there will be many resonances within a loop and 
smaller length scales will be created due to phase mixing, causing a cascade of energy to 
smaller scales where dissipation becomes more efficient. 



In the follow up paper by 



Tomczyk &: Mclntoshl (120091 ) more detailed analysis of the 



same CoMP data set showed th at propagating way es tra yeling a larger distan ce suffered 



greater damping. This inspired 



Pascoe et al. 



torn and 



Terradas et al. 



2010l ) to model 



propagating ki nk waves damp e d by the process of resonant absorption. Most relevant to the 



present study. 



Terradas et al. 



(J2010l ) derived simple analytical expressions for the damping 
length of kink waves, showing damping length is a monotonically decreasing function of 
frequency, the TGV relation. This has a consequence that solar waveguides with transverse 
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inhomogeneity, e.g., coronal loops, will act as low pass filters for propagating kink waves. It 



is the main purpose of t 



by 



Tomczyk et al. 



lis Letter to investigate if the propagating Alfvenic waves observed 



(120071 ). realistically interpreted as the kink mode from MHD wave theory 



exhibit this frequency dependent damping. 



Theoretical model 



Traditionally, since standing kink waves in coronal loops were ob served using TRAC E 



the m easu rement of daniping t imes was of primary interest, see e.g.. 



(l2003|)and 



Aschwanden et al. 



Arregui et al.l (120071 ). However, for propagating kink waves it is more appropriate 



to study the damping length Ld, with the expression for wave amplitude given by 



A{s) = Aoexp{ —— 



where s is the distance along the waveguide and Aq is initial amplitude. From 



(2) 



Tomczyk et al. 



(120071 ). the transverse scale of propagating disturbances has an upper limit of about 9 Mm 
and wavelengths are A ^ 180 Mm, so assuming that the average waveguide observed with 
CoMP is a fluxtube of radius R ^ 4.5 Mm, ki nk waves are in the lo ng wavelength regime 
where R/X ^ 1. In this limit it was shown by 
is simply given by 



Terradas et al. 



mm that the TGV relation 



Lb — fph'^ET, 



(3) 



where / is the frequency defined by / = Vpi^/X and ^e is an equilibrium parameter 
dependent on length scale of the density inhomogeneity. Equation ([3]) demonstrates that 
Lj) is inversely proportional to /, i.e., the rate of damping per unit length will be greater 
for higher frequency waves than low frequency waves. The parameter ^^ can be calculated 
precisely for a chosen equilibrium model, e.g., if we choose a thin inhomogeneous boundary 
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layer with a continuous sinusoidal profile decreasing between pi and pe then 

2Rpi + Pe 

IT I Pi- Pe 

where / is the thickness of the boundary l ayer (see e.g. , 



(4) 



Goossens et al. 



1992 



Ruderman fc Roberts 



2002 : 



Goossens et al 



2OO2I ) . Hence Equations ([3]) and 



demonstrate that the efficiency of damping due to resonant absorption also depends on the 
thickness of the non-uniform layer and the steepness of the transverse gradient in density. 



From Figure 1 (Panel A) in iTomczyk fc Mcintosh! (120091 ). it can be seen that the 



Tomczyk fc Mcintosh 



i ntegr ated wave path in the coronal loop system is off limb, with 
(I2OO9I ) assuming the path is in the plane of sky. The approximate semi-circular annular 
geometry of the this wave path is shown in Figure (H with the direction of outward and 
inward propagating waves denoted by arrows. We show the analyzed wave paths by the 
shaded grey region, whic h consists of many wayes tra veling along different coronal loop 



structures, integrated by 



Ratio fSNR), c.f.. Fig ure 4 



(Panels B and C) in 



Tomczyk fc Mclntoshl ( 20091) to improve the Signa 



Tomczvk fc Mclntoshl (120091 ). Although 



to Noise 



Tomczyk fc Mclntoshl (120091 ) did not integrate the entire loop lengths, i.e., from the outward 
footpoint (s = 0) to the inward footpoint {s = 2L), the important point is that in the 
grey region, inward waves have traveled further and therefore the damping will be greater 
relative to the outgoing waves. 

Let the initial outward at s = and inward power at s = 2L for frequency / be denoted 
by Poutif) and -Pin(/) respectively. Then since P{f) oc v4^(/), by Equations ([2]) and (jH]) the 
spatially averaged outward power in the shaded grey region s E [0,L] for frequency / is 

2/ 



^ph6 



and for the inward power is 



s ds 



^ph6 



s I ds. 



(5) 



(6) 



Note that for Equation ([6]) integrating the power in the grey region from s = L to s = for 
the inward wave is equivalent to integrating from s = L to s = 2L for the outward wave. 
Defining (P(/)),^tio = (^(/))out / (^(/))in' i^ can be shown by Equations (P and © that 



(^(/))ratic 



Poutif) 



exp 



2L 



-/ 



In the next Section we shall do an expone ntial least squares fit using E quation ([7]) with 



(7) 



CoMP data from the same observation by 



Tomczyk fc Mclntoshl ([20091) in the frequency 



range — 4 mHz where the SNR is strongest. 



4. Least squares fit to data 

In Figure [2] it can be seen that the CoMP data (see triangles) shows a clear trend of 
(-P(/))ratio increasing with frequency /, i.e., that higher frequency waves are damped more 
than their low frequency counterparts. We add the caveat that in the CoMP data, this 
trend continues to about 8 mHz, then there is a turn over in the ratio {P{f)) j-atio- -^^ ^^^ 
present time it is unclear if this is simply due to background noise or is the result of some 
other physical mechanism/s. Understanding this high frequency trend should be the focus 
of a future study. However, there is a strong level of confidence in the trend of (-P(/))i.atio 
in the range — 4 mHz, which can be explained by the theory of propagating kink wave 
resonant damping. To illustrate this we implement a least squares fit of the power ratio 
function given in Equation ([7D to the CoMP data in this range. From the estimates of 
Tomczyk fc Mclntoshl (120091 ) our fixed parameters are fph = 0.6 Mm s~^ and path length 
L = 250 Mm for both outward and inward waves, leaving the free parameters as Pout /-Pin 
and ^E- It is found that the best estimates are Pout/-Pin = 0.91 and ^e = 2.69 (see solid 
line in Figure 2) with 95% confidence bounds Pnnt/Pm = 0.67 — 1.8 9 and ^e = 1.15 — 3.49. 



Although the wave paths analysed by 



Tomczyk fc Mclntoshl ( 120091 ) are integrated over 



many coronal loop structures, the spatially averaged estimate of .^e is consistent with 
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outward 
direction 



s = L 



inward 



5=2/. 




5 = 



Fig. 1. — Simple cartoon of the observed semi-circ ular annulus geometry o 



the coronal loop 



system. The integrated wave paths analyzed by iTomczyk &: Mclntoshl ( 120091 ) (shaded in 
grey) are approximately only half the length along the total loop system [L ^ 250 Mm) 
since the data in this segment had the best SNR. The direction of outward and inward wave 
propagation is shown by the solid and dashed lines respectively. 




Fig. 2. — Power ratio against frequency best fit using Equation ([7]) with CoMP data. Fixed 
parameters are f ph = 0.6 Mm s~^ and L = 250 Mm. Best estimate parameters are Pout/ Pin = 
0.91 and .^e = 2.69 (solid line). The 95% confidence intervals for the simultaneous functional 
bounds are shown by the dashed lines. Note that in the CoMP data SNR decreases with 
increasing frequency. 
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the dampe d standing; kink waves in individual coronal loops observed by TRACE where 



^E ~ 1 — 4 (JAschwanden et al 



20031 ). It is interesting to note that the average half length of 
loops estimated with TRACE data is about L = 110 Mm due to its restricted field of view. 
The longer value of L = 250 Mm observed with CoMP suggests that the rate of damping 
is independent of loop length, i.e., long loops have similar transversal length scale density 
inhomogeneities as sho rt loops. This has also been independently confirmed recently by 



Verwichte et al. 



(|2010[ ) in the study of a long loop of L = 345 Mm using combined SoHO, 
TRACE and STEREO observations, where they found ^e = 1-48 for a post-CME standing 
kink wave. Thus, the analysis of observational cases so far, suggest they are all within the 
valid regime of the thin tube approximation of resonantly damped kink waves. 

Another interesting feature shown in Figure [2] is the estimated least squares value of 
Pout /-Pin ~ 1, suggesting that power generated at both the outward and inward footpoints is 
approximately equal. The 95% confidence intervals for the simultaneous functional bounds, 
i.e., calculated with all predictor values, are also shown in Figure [2] by the dashed lines. It 
can be seen that there is a trend of decreasing confidence for / ^ 3 mHz, most likely due to 
the fact that the SNR for the CoMP data decreases with increasing frequency. 

From estimates of L, Vp^ and ^e, using Equations ([2]) and ([3]) we can calculate the 
percentage power loss as a function of frequency for the observed propagating kink waves 
after one travel time between the two loop footpoints, s = and s = 2L (see Figure [3]). The 
kink waves with / ^ 2.5 mHz lose at least 50% of their power and as explained previously 
in Section [21 this broadband frequency power is converted into Alfvenic azimuthal motions 
in many resonant surfaces. 
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Conclusions 



In this Letter we establislied more evidence for the damping mechanism of resonant 
absorption by analyzing observational data from the CoMP. Crucially, this instrument has 
established in the solar corona there are ubiquitous propagating low amplitude (^1 km s~^) 
Alfvenic waves with a wide range of frequencies. Realistically interpreting these propagating 
waves as kink modes, it was predicted they should exhibit a frequency dependent damping 
length due to resonant absorption, governed by the TGV relation which shows that 
naturally occurring transversal plasma inhomogeneity in coronal magnetic flux tubes causes 
them to act as low-pass filters. It was found the observed frequency dependence on damping 
length can be explained by the kink wave interpretation, at least up to about / = 8 mHz, 
and furthermore, the spatially averaged equilibrium parameter describing the length scale 
of transverse plasma density inhomogeneity over a system of coronal loops is consistent 
with the range of values estimated from TRACE observations of standing kink modes. 
Due to its restricted field of view, full loops observed with TRACE are shorter than the 
average loop length of the coronal loop system observed with CoMP. The similar inferred 
transversal plasma inhomogeneity length scales, suggest this parameter is independent of 
loop length. Importantly, it was also found that the estimated least squares value of the 
ratio of initial outward and inward power was approximately unity, suggesting the kink 
wave power generated at both footpoints in the lower atmosphere was almost equal. 

Previously much work has been undertaken to study the properties of the well known 
standing kink modes in post-flare coronal loops. The results of this Letter now open the 
way to investigating the wave contribution to coronal heating by exploiting the latest high 
spatial/temporal resolution observations of propagating kink waves modes. An accurate 
estimate of coronal kink wave damping rates due to resonant absorption is crucial in 
order to quantify the amount of this wave energy that can contribute to plasma heating. 
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Broadband frequency disturbances propagating along transversally inhomogeneous coronal 
flux tubes will generate many resonances that will phase mix, causing a cascade of energy 
to smaller scales where dissipation becomes more efficient, thereby creating the necessary 
conditions for plasma heating. Unlike the larger amplitude standing kink mode coronal 
loop oscillations which need relatively high energy but rare excitation events like flares 
or CMEs, small amplitude propagating kink waves are ubiquitous. This makes them an 
ideal candidate for study with coronal observing instruments of sufficient Doppler velocity 
resolution such as CoMP, Hinode and the planned Advanced Technology Solar Telescope 
(ATST). 

The authors thank S. Tomczyk for his help and expertise regarding the CoMP data. 
J.T. acknowledges the Universitat de les Illes Balears for a postdoctoral position and 
the funding provided under projects AYA2006-07637 (Spanish Ministerio de Educacion y 
Ciencia). M.G. and G.V. acknowledge support from K.U.Leuven via GOA/2009-009. M.G. 
also expresses gratitude to both the FWO Vlaanderen and UIB. 
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Fig. 3. — Percentage power loss against frequency of a kink wave after one travel time 
between the two loop footpoints s = and s = 2L, a distance of 500 Mm. The dashed lines 
indicate the percentage power loss up to the 95% confidence bounds of estimated parameter 
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